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We present the first field-theoretic calculations of the contribution made by cosmic strings to the 
temperature power spectrum of the cosmic microwave background (CMB). Unlike previous work, 
in which strings were modeled as idealized one-dimensional objects, we evolve the simplest example 
of an underlying field theory containing local U(l) strings, the Abelian Higgs model. Limitations 
imposed by finite computational volumes are overcome using the scaling property of string networks 
and a further extrapolation related to the lessening of the string width in comoving coordinates. 
The strings and their decay products, which are automatically included in the field theory approach, 
source metric perturbations via their energy-momentum tensor, the unequal-time correlation func- 
tions of which are used as input into the CMB calculation phase. These calculations involve the use 
of a modified version of CMBEASY, with results provided over the full range of relevant scales. We 
find that the string tension fi required to normalize to the WMAP 3-year data at multipole i = 10 
is Gjj, — [2.04 ± 0.06(stat.) ± 0.12(sys.)] x 10~ 6 , where we have quoted statistical and systematic 
errors separately, and G is Newton's constant. This is a factor 2-3 higher than values in current 
circulation. 



I. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) radiation have helped establish a strong case for 
cosmic structure to have grown from primordial pertur- 
bations created via inflation. On the other hand, theories 
in which the seeding of structure is primarily attributed 
to the presence of topological defects [lj have failed to 
match the data. These defects, of which cosmic strings @ 
are the prime example, have therefore been relegated to 
at most secondary phenomena with CMB measurements 
providing upper limits on their relative importance. Al- 
though the precise constraints given by current data de- 
pend upon the defect model and the data sets chosen, 
recent calculations have shown that their maximum al- 
lowed effect upon the CMB temperature anisotropics is 
at a level of around 10% in the power spectrum Hi- 

Nevertheless there is once again great interest in topo- 
logical defects, particularly in the case of cosmic strings, 
thanks to recent theoretical and observational results (see 
Q for a review). Cosmic strings seem to be viable post- 
inflation entities in supersymmetric grand unified theo- 
ries (GUTs) 0, superstringtheory @,©, EE El El and 
hybrid inflation scenarios [13L EH- They have also fea- 
tured as a possible explanation for the lens candidate 
CSL-1 [IE EE E3, El], which initially appeared to be 
two images of the same galaxy, without the arc-like dis- 
tortions that would result from lensing by a spheroidal 
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distribution of matter. While new HST data [IE] has 
revealed that CSL-1 is actually two interacting galaxies, 
this story now serves to highlight one means of cosmic 
string detection. In a second observational case, an os- 
cillating loop of cosmic string has been discussed as a 
possible explanation for the synchronous brightness fluc- 
tuations in the two images of the (normal) gravitational 
lens system Q0957+561 HO]. 

In the present work we consider the effects of an en- 
tire network of cosmic strings on the CMB, with the hope 
that they may be detected by future measurements. Such 
observations would allow inferences to be made about 
the network properties, which in principle depend on the 
nature of the underlying theory and so provide an im- 
portant window on physics at very high energy. Notably, 
we present the first CMB calculations for cosmic strings 
to employ simulations of a field theory, in this case the 
Abelian Higgs model, which contains (local) U(l) strings. 
This is a considerable computational challenge that has 
been made possible via the use of two extrapolations, 
each justified by the results from the simulations them- 
selves. These enable us to ascertain the contribution 
to the CMB power spectrum from (traditional) cosmic 
strings and to compare our results with previous deter- 
minations, discussed momentarily, which have involved a 
greater degree of modelling but have the advantage that 
they are less computationally demanding. 

The strings are considered to have formed at the end 
of inflation or in a later phase transition. They would 
then additionally perturb the cosmic fluids at all sub- 
sequent times and there would hence be two sources of 
anisotropy in the CMB. These may in fact be taken to 
give independent contributions to the CMB power spec- 
trum, and therefore here we calculate that from strings 
which is to be added to the dominant contribution from 
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inflation. Even if the strings are directly related to the 
inflation mechanism, their complex non-linear evolution 
destroys correlations with earlier times. Further, the ho- 
mogeneity of the CMB implies that the effect of the in- 
flationary perturbations on the strings will be negligible 
and can be ignored in linear cosmological perturbation 
theory. Therefore, the perturbations in a quantity X(k) 
arising from the two mechanisms can be calculated sep- 
arately and the power spectrum given by: 



x*x) = (X^jti 



XtefXdei) + 2 (X* nf Xdet) , ( 1 ) 



but with the cross-term simply zero. The normalization 
of the string component is then a free parameter related 
to the energy-scale of the model or equivalently the ten- 
sion /x, normally expressed as the dimensionless combi- 
nation Gfi. While we present our calculation method 
and results for the CMB power spectrum here, the pre- 
cise constraints on G\x will be addressed in a separate 
publication [21 1. 



A. Previous cosmic string power spectra 

All previous determinations of the CMB power spec- 
trum for local cosmic strings have relied upon modelling 
the strings as idealized objects of infinitesimal width. In 
reality, the string width is related inversely to the energy 
scale present in the corresponding theory, which could be 
as high as the grand unification scale 10 16 GeV, whereas 
the string separation is comparable to the Hubble radius 
(corresponding to 1CP 42 GeV in the present epoch). The 
idealization is hence justified by the fact that, at times of 
importance for CMB calculations, the difference in scale 
is enormous. 

Past results have stemmed from employing this ap- 
proximation using either (i) Nambu-Goto simulations of 
connected string segments [22|, [H, HH HE] \ or (ii) a model 
involving a stochastic ensemble of unconnected segments 
0; US H3j Hl| ■ The first case involves simulations of the 
dynamical equations for ideal relativistic strings, either 
in a Friedman-Robertson- Walker (FRW) [H, [H[ or 
Minkowski [24| space-time. In (ii), the segments are ran- 
domly selected for removal so as to give sub-horizon de- 
cay, with the particular parameters of the model chosen 
to provide a match to, for example, the segment density 
seen in simulations. Although this approach involves a 
greater degree of modelling, the CMB results match the 
form given by the more computationally intensive sim- 
ulations of Contaldi et al. [24]]. However, it should be 
pointed out that the Contaldi et al. [24[ (and derived 
papers [29|, [3(|) present the only previously published 
simulation-based CMB calculations for the very angular 
scales (£ — 300 — 700) at which the string contribution is 
shown in to peak in that reference. 

Unfortunately, there have been questions raised as to 
the accuracy of the Nambu-Goto simulations themselves 
with regard to small-scale structure, loop production and 



the decay of strings. For example, in the Minkowski 
space-time simulations using the Smith- Vilenkin algo- 
rithm (as used by Contaldi et al. [13]), the only means 
via which the network can lose energy is in the removal 
of small loops, created by self-intersection events. In the 
conventional picture, these loops would decay into grav- 
itational radiation. However in those simulations, loop 
production is seen to occur most frequently at the small- 
est scales involved {31], which is the lattice grid in the 
Smith- Vilenkin algorithm. Recent simulations using dif- 
ferent algorithms in Minkowski space (321 . [33[ or in an 
FRW universe [H, [35| point to the initial correlation 
length as the dominant scale of loop production, although 
one group claims some sign of evolution toward a scaling 
form for the loop production function [35| . 

It may be the case that the true Nambu-Goto physics 
involves increasing energy density in loops of smaller and 
smaller size until the Nambu-Goto approximation fails 
close to the string width. If so, then particle production 
may be more important than gravitational radiation as 
the means of string decay. A second possibility is that 
loop production peaks at some small fraction of the hori- 
zon scale and it is simply that this has yet to be observed 
in even the largest Nambu-Goto simulations to date. 



B. Power spectra from field simulations 

In contrast to the above, the present approach to CMB 
calculations for cosmic strings employs simulations in 
which the string width itself is resolved. That is, the 
Abelian Higgs model (in the classical approximation) is 
represented on a lattice and the strings are evolved in 
terms of their constituent fields. Although this is the 
simplest model to exhibit local strings, available compu- 
tational resources limit the size of the simulation volume 
to a few hundred times the string width. Given the above 
horizon size considerations, the epoch of interest for the 
CMB cannot be simulated directly. 

Fortunately, cosmic strings are believed to evolve to- 
ward a scaling regime, in which the network appears sta- 
tistically the same at all times relative to the causal hori- 
zon. This enables the statistical results from a small re- 
gion to be scaled up to larger volumes at later times, 
which is precisely that required for CMB calculations. 
However, it is breached at the radiation-matter transi- 
tion and hence it is the case that scaling solutions in 
both eras must be studied. 

In fact, scaling has previously been employed in all 
high resolution power spectrum calculations (I > 300), 
either in unconnected segment models to control string 
decay or to boost the range of scales over- which Nambu- 
Goto string simulations can provide data. Such ideal 
simulations are scale-free, however the incorporation of 
the string width in the present case means that scaling 
is used to explicitly translate the results to larger scales 
and as well as boosting the dynamic range. 

Studies of the Nambu-Goto dynamics and of the field 
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FIG. 1: Slices through a 512 simulation in the radiation era, showing the Abelian Higgs analogue of magnetic flux density. 
Magnetic flux tubes run along the cosmic strings, which appear as dark regions, with a shape dependent upon the nature of 
the string intersection with the slice. Varying left to right, the horizon size (measured as 2r) is 0.35, 0.63 and 1.00 times the 
box-side while the string width decays inversely with the scale factor a and hence as r^ 1 (a a t). Note the decay products 
visible in these images. 



micro-physics provide complementary illuminations of 
cosmic string dynamics and while field simulations have a 
smaller dynamic range, they have a number of attractive 
properties, particularly with regard to CMB calculations. 
Firstly, the decay products of the strings give rise to addi- 
tional CMB perturbations but these must be included in 
an ad hoc manner in Nambu-Goto simulations. Energy 
is transferred via intersection events from long strings 
down to small loops, which are traditionally considered 
to decay into gravitational radiation and are therefore re- 
moved from the simulations. It is not the case that grav- 
itational waves are likely to be included self-consistently 
in any simulation in the near future, due to both the nu- 
merical complexity, and technical problems arising from 
computing the back-reaction of the gravitational radia- 
tion on the string network. Contaldi et al. [2J| therefore 
included a compensating fluid via which energy was con- 
served. On the other hand, particle production is mod- 
eled automatically in classical field simulations in terms 
of oscillations in the fields, which can be seen in Fig. [TJ 

Secondly, field simulations allow a wide class of mod- 
els to be studied which cannot be approached via the 
Nambu-Goto approximation. For example, we have al- 
ready employed the method described here to the ad- 
ditional case of semi-local strings [36(, with the present 
work enabling a comparison between cases without ques- 
tion marks over possible systematics concerning the two 
simulation approaches. 

Thirdly, a classical field theory simulation is a well- 
controlled approximation: the H — > limit of the un- 
derlying quantum field theory, and given enough com- 
putation even the quantum corrections could be com- 
puted using stochastic quantization techniques [37] . The 
Nambu-Goto approximation is itself an approximation of 
the classical field theory, and describes only smooth field 
configurations representing strings with large radius of 



curvature relative to their width. It fails at cusps and 
kinks [H, Ht| which are inevitable features of a mobile 
and intersecting string network. 

Field simulations have previously been involved in 
CMB calculations for the case of global defects [53, [5H . 
In the global case, the defect cores contribute a small 
fraction of the energy compared to the large-scale field 
variations and therefore the cores may be left unresolved, 
with the non-linear er-model used to approximate the dy- 
namics. However, in the local string case, the string cores 
are an important source of energy and momentum and 
they must be resolved on the computational grid. This 
presents a serious problem: the string width is a fixed 
physical length scale, which rapidly decreases in the co- 
moving coordinates used to represent a region of an ex- 
panding universe. Either the strings are too narrow to 
be resolved at the end of the simulation, or so large ini- 
tially that a representative network cannot be formed. 
This is illustrated in Fig. [TJ where it is apparent from 
the left-most image that, in this large simulation during 
the radiation era, the string width is not greatly smaller 
than the string separation. Strings cannot form greatly 
earlier but the horizon is already around one-third of the 
box-side. And that is despite the string width being close 
to the resolution limit at the end of the simulation. In the 
radiation era the cosmic scale factor a is proportional to 
conformal time t, but in the matter era it varies quadrat- 
ically: aocr 2 . This makes the effect even more of a prob- 
lem in the matter era and the available range in t is very 
limited. 

It is hence required to adopt an approach similar to 
that used in the Abelian Higgs simulations of Moore et 
al. [42j (and domain wall simulations of [43]), in which 
the equations of motion were artificially modified such 
that the string (or domain wall) width became comov- 
ing. Focusing on the Abelian Higgs case, Moore et al. 
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have performed simulations in the radiation era, both 
with and without this artificiality, and report little dif- 
ference in the dynamics. Unfortunately, the approach 
introduces a breach of the conservation law for the very 
energy-momentum tensor through which the strings in- 
teract gravitationally with the cosmic fluid. While not 
especially relevant for the work of Moore et al., this is 
potentially important for CMB calculations. 

The approach used here is therefore to modify the 
equations so that the comoving string radius r varies with 
the scale factor as: 



(2) 



such that s controls the sensitivity of the string width 
upon a. We then make CMB calculations using values of 
s between s = (comoving width) and the closest value 
to the true case (s = 1) that our facilities permit. Hence 
in principle, the effect of s upon the CMB results can be 
ascertained and accounted for. 

CMB predictions are made using the unequal-time cor- 
relator (UETC) approach [4fJ, |4J|, following the Durrer 
et al. [4l|, |45| formalism. In this method, the statistical 
information taken from the simulations are the two point 
correlation functions of the energy-momentum tensor: 



U\KH„(k,T, t') 



(3) 



between unequal times r and t' . Although this approach 
limits the calculation of CMB correlation functions to 
power spectra, it enables a modified version of an infla- 
tion CMB code, CMBEASY @ in the present case, to 
calculate the perturbations created by the defects. As 
will be explained in Sec. IIII1 this is via a re-expression 
of the UETCs as a sum of coherent source functions 
which drive the string perturbations. Scaling, statisti- 
cal isotropy, and causality constrain the form of UETCs 
and their dependence upon the absolute times is trivial. 
It is merely the ratio r/r' that is important and this en- 
ables their application over a large range of times, and 
therefore CMB scales. This is in contrast to alternative 
approaches [HI, HE H3] , which afford the production of ac- 
tual CMB maps, useful for non-Gaussianity studies, but 
are limited to only a small range of scales due to compu- 
tational constraints. However, the UETC approach still 
requires the simulations to provide data for the range 
of time ratios that the UETCs are non-negligible. While 
UETCs decay for large or small ratios, the range of times 
over which the strings can be represented is increasingly 
limited as s is increased. Since the strings take some 
time after formation to reach scaling to the accuracy re- 
quired, it is this range of time ratios required that limits 
the practical value of s. 



II. STRING SIMULATIONS 

A. String model 

The Abelian Higgs model, involving a complex scalar 
field 4> and a gauge field A^, has Lagrangian density: 



C 



4e 2 M 



A 



(4) 



Here = d^+iA^ is the gauge-covariant derivative and 
Ffj, v = d^A u — d v A^ is the field strength tensor, with e 
and A dimensionless coupling constants. Note that rela- 
tive to many references, the gauge field is rescaled using 
e, which proves useful when controlling the dependence 
of the string width upon a. 

The model obeys the local U(l) symmetry: 



(j) — > <j) exp(za) 
A^l ► A^ d^a, 



(5) 
(6) 



which is broken in the low-temperature regime (without 
quantum corrections to the potential) by any choice of 
vacuum \ip\ — (f>Q. That the vacuum manifold is a closed 
loop allows for topologically stable string defects, around 
which the phase of cf> has a net winding of 2irn {n £ Z ). 
Although the phase is not a gauge-invariant quantity, 
such a winding would require an infinite gauge transform 
for its removal. However, the gauge field counteracts the 
gradient energy associated with the winding such that, in 
contrast to global strings, the energy of the configuration 
is localized to a comoving radius: 



1 



iVA (j)Q ' 



(7) 



within which \<fi\ approaches zero. As the gauge field at- 
tempts to counter the phase gradients near the core, it it- 
self acquires a significant curl around the string, resulting 
in a magnetic flux tube which traces the string. For the 
present work we adopt the Bogomol'nyi ratio A/2e 2 = 1, 
in which the characteristic scales of the magnetic and 
scalar energies are equal. It is important to note that de- 
spite the analogy with electromagnetism, in the broken- 
symmetry phase, perturbations about the vacuum reveal: 

(8) 



and that the gauge field is massive. This model contains 
no massless modes and therefore all group velocities are 
less than the speed of light, except for waves along the 
strings themselves. Hence, the normal causal run-time 
limit imposed by the simulation boundaries may be ex- 
tended as the strings are highly curved and therefore no 
disturbances can traverse the box by this time. 

Variation of the action corresponding to Eq. (U|) and 
making the gauge choice Aq = yields the dynamical 
equations for a flat FRW space-time as: 



,A 



■-oj 



f) F- 



-2a 2 e 2 lm\<b*Di 



(9) 
(10) 
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The notation here is such that over-dots denote differ- 
entiation with respect to conformal time r, and di with 
respect to comoving Cartesian coordinates. Further, the 
system obeys the constraint: 

-d l F 0l = -2a 2 e 2 lm[ ( j ) *<Pl (11) 

which is analogous to Gauss' law, and thereby conserves 
the 4-current density: —2a 2 e lm[(f)* D^ifi). 

It is the presence of the a 2 factors on the right-hand- 
sides of these equations that causes the string width to 
lessen in comoving coordinates and causes the problem 
highlighted in the previous section. Moore et al. 
have evolved Eq. © and (fTT)]) in the radiation era and 
then additionally with the factors of a 2 removed from 
the right-hand-sides. The result was that not only did the 
string length show scaling but that the scaling forms were 
very similar, allowing them to use the modified equations 
to study the matter era also. However, it is the desire 
here to have a controllable dependence of the width upon 
a, parametrized by the parameter s according to Eq. ^ . 
One means of achieving this is to raise the dimensional 
coupling constants to time-dependent variables according 
to: 

A = Aoa~ 2(1 - s) (12) 
e = e a-^- s \ (13) 

Variation of the action now yields the artificial equations 
of motion: 

4> + 2-j - DjDjtfi = -a 2s ^ (H 2 - 2 )0 (14) 

F oj + 2(l-s)-F 0j - diFij = -2a 2s el lm[<t>* Drf] , (15) 

and the constraint: 

-9 lJ F Ol = -2a 2s e 2 Xr72[0*0]. (16) 

Notice that there is an additional damping term for s < 1 
that is created by the lowering of e with time, and is 
required for the dynamical equations to preserve Gauss' 
law in this consistent form. Therefore this approach, even 
in the case of s = 0, is not precisely the same as that of 
Moore et al. [13]. 



method of Moriarty et al. [49j . Rather than discretizing 
the dynamical equations directly, which does not gener- 
ally lead to an algorithm that preserves the Gauss con- 
straint, the Moriarty et al. approach is to discretize the 
Hamiltonian. Then discrete equations of motion, includ- 
ing a discrete version of the constraint equation, can 
be derived from it. This is broadly the approach fol- 
lowed here, except that the action is discretized rather 
than the Hamiltonian, preferable given that these simula- 
tions are performed for an FRW rather than a Minkowski 
metric, and additionally the variables A and e are time- 
discrctized. 

The discretization preserves the gauge symmetry in the 
form: 

^x,r ^ exp (j a x,T) 

A^ +h - - ± (a x - +1 - a^) (18) 

Af ilT -» Af iS ' T - ^- (a x+ l T - a x ' T ) . (19) 

The notation is such that A* + ^' T signifies that the spa- 
tial components of the gauge field are represented on the 
links half-way between lattice steps. This creates a con- 
sistent centered-derivative transform that is in-line with 
the geometric interpretation of the gauge field: that it 
performs a local rotation of the field-coordinates to form 
the gauge-covariant derivative. Hence, this derivative is: 

(D^T^ = J- (>+J - expH^V) , (20) 

with: 

(9* + * j = AxA* +ii , (21) 

making the gauge-field a purely angular variable. 
The entity: 

A x +i i+ij = ^ i+ i_ d ^_^ i+i _ e ^ (22) 



B. Lattice discretization 

It has become standard in the literature [13, Ell to 
discretize the Abelian Higgs model on a lattice via the 



is invariant under the discrete gauge transform and is 
used to build the Fy magnetic term in the action. How- 
ever, following Wilson [501 ] . this is performed in a manner 
that preserves this angular nature and the total action 
integral is represented as the discrete sum: 



s = At (^EE (- s^tEE [' - -OF*")] + 2 ^ (A L, w E - 
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(I0 ; 



X,T|2 



Hence the first term tends to ^FijFij in the continuum 
limit since 1 — cos a; — > ^x 2 as x — > 0. 

In the temporal gauge Aq = used for evolution, there 
is no requirement to treat the Foi electric term in such a 
manner and the second term is essentially just the time 
derivative of the gauge field. Note however that this term 
involves a centered and therefore mid-step derivative and 
so e has likewise been referenced mid-step. The situa- 
tion is then very similar for the third term, which (in 
this gauge) is just the time derivative of the scalar field, 
with the scale-factor referenced mid-step. The remaining 
terms then involve the spatial gauge-covariant derivative 
and the potential. 

Variation of this sum with respect to X ' T or (9 x+5l ' T ; 
for a given x and t, then yields the evolution equa- 
tions. The constraint may be derived by including Aq 
(although there is no requirement to do so via the full 
Wilson method) and then minimizing the action with re- 
spect to it. Preservation is assured (to machine precision) 
in the temporal gauge due to the discrete gauge symme- 
try and that the constraint results from an independent 
minimization of the action sum. The dynamical equa- 
tions in the chosen gauge allow for a leap frog update 
algorithm, in which the fields are updated using their 
stored time derivatives, and then the time derivatives 
at a site are updated using the field values at that site 
and its neighbours. There are no off-site references made 
to the quantity being updated and therefore there is no 
additional temporary storage required. The evolution 
hence involves 10 single-precision floating point numbers 
per site, or for an N 3 lattice, 40 (TV/1024) 3 GB of data 
updated each timestep. 

The resolution of strings at the end of a simulation is 
assured by setting a en d = 1, Ao = 2, eo = 1, and then 
Ax = O.50Q 1 , which has been shown to yield good string 
resolution in Minkowski space-time simulations (a = 1) 
. The choice of At /Ax — 0.2 gives good covariant 
energy-momentum conservation in the s — 1 case (in the 
radiation era). For example during an s — 1 phase in 
512 3 simulations, the measured decrease in the mean en- 
ergy density was found to be within ±7% of that expected 
from Hubble damping. For s = 0.3 in the radiation era, 
the rate of energy loss is 10 — 20% greater than the con- 
servation law would predict. As discussed in the pre- 
vious sub-section, the absence of massless modes allows 
the simulation to be run longer than the half-box cross- 
ing time for light, with T en( j = 1.25(iV/2)Ax used here 
(for s < 1 cases). However, these late times are used to 
provide the least important UETC data. A single sim- 
ulation requires ~ 500 CPU-hours of processor time on 
the UK National Cosmology Supercomputer [5l|. Paral- 



lel processing is afforded by MPI via a library, written 
as part of this work, which allows the rapid development 
of parallel field-evolution simulations using a convenient 
C++ object-orientated interface (52j. 



C. Initial conditions 

An initial field configuration consisting of a scaling net- 
work of strings is difficult to achieve directly. Previous in- 
dependent simulations of Abelian Higgs string networks 
[42I HH have used different initial conditions and have 
recovered consistent results for the string length density 
in Minkowski space-time during the scaling regime. Ac- 
cepting that the precise nature of the initial conditions 
is unimportant, it is desirable that they give a rapid con- 
vergence toward scaling but with minimal complexity in 
their generation. 

The initial field configuration must satisfy the discrete 
form of Gauss' law and the simplest means of achieving 
this is to set all temporal derivatives and the gauge field 
to zero across the entire lattice. Then <j> is set to have 
modulus <pQ with a random phase assigned to each lattice 
site. With the initial time set such that Ax is of order 
the causal horizon, this sets up independent phases in 
each initial horizon volume and is a rough approximation 
to the results of a phase transition. The system is then 
evolved according to the discrete equations of motion and 
Hubble damping relaxes the system into a network of 
cosmic strings. This approach is in contrast with Moore 
et al. 42] in which the network relaxation is achieved 
using a period of diffusive (first-order) evolution, with the 
second derivatives removed from the equations of motion. 

However, strings cannot form until the causal horizon 
is larger than their characteristic width, since otherwise 
gradients exist on scales so small that there is sufficient 
energy for the entire volume to rise up the potential to- 
ward |0| = 0. Although the choice of s < 1 lowers the 
initial string width relative to its final value, strings will 
still not form for some time. Since the network will take 
further time to reach the scaling regime, it is desirable 
to accelerate the formation process so that UETC data 
may be taken over as large a ratio of times of possible. 
Hence the parameter s is in fact set to be negative prior 
to a time r s such that the string width in fact grows from 
an initially low value, before s takes on its final positive 
value and the width shrinks during the second phase of 
the simulation. The string width hence varies according 
to Fig. [2] and a network of strings forms by r ~ lO^" 1 . 
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FIG. 2: The variation of the string width r and coupling 
parameters with conformal time r (in units of (fig 1 ) for the 
s = 0.3 simulations described in Section Til CI The subscript 
indicates the value at the end of the simulation. 



III. CMB CALCULATION METHOD 
A. UETC approach 

The evolution of the cosmological perturbations may 
be described by a linear differential equation of the fol- 
lowing form: 

-D a c(k,a,a,p,...)X a {k,T') = S c (k,T'), (24) 

with X a denoting the Fourier Transform of X a . The 
linear differential operator T> a b includes quantities from 
the background FRW universe such as the mean physi- 
cal density p and the cosmic scale factor. It acts upon 
the metric, matter and photon perturbations described 
by the vector X a , with the source term S c describing 
the active seeding due to the defect presence — their 
energy-momentum tensor. For the homogeneous case 
corresponding to inflation (S c = 0), this equation set 
can be solved using the standard CMBEASY code. In 
principle therefore, if 5 c (k, r') is known, then this inho- 
mogeneous set can be solved using a Green's function 
Qac(k,T,T r ) to give the perturbation power spectra for a 
wave vector k and time r as: 

(x a Xt) = JJ dr'dr"g ac {r')Gt d (r") {s c {r)S* d (r")) . 

(25) 

Although this is not the actual method used here, this 
equation shows that the data needed from the simulations 
for CMB power spectra calculations are the bracketted 
terms on the right. These are the Fourier transforms of 
the two-point correlation functions: 

EUy,T,r') = i / d 3 x <S c (x,r)S d (x- y ,T')), (26) 



with the normalization of the Fourier transform chosen 
as: 

U cd (^ r, r') = ^J d 3 y U cd (y, r, r') e~^. (27) 

Here V is the fiducial simulation volume and its inclusion 
yields U cd with the same dimensions as S c (and so T^ v ) 
squared: [U cd ] = [U cd ] = (time) -4 . 

For scaling sources, a statistical measure of the dynam- 
ics should be dependent upon a single scale d. While the 
energy-scale 4>o should not affect the spatial distribution 
of strings, it does set the normalization of the energy- 
momentum tensor as <pQ. Hence assuming scaling, C4 C 
may be written as: 

= 5-7)' (28) 

with / a dimensionless function. The scale d must be 
symmetric in the two times involved and hence may be 
written as: 

d = VTr 1 9 (t/t') = \ZVT' g(r' It). (29) 

In this form the final two inputs in / provide the same 
information, and further the dimensionless function g can 
be absorbed to yield: 

= ^ F (j=,Ly (30) 

The Fourier transform then gives: 

U cd (k,r,r) = rtlc(kV7^,T/T>). (31) 
yrr V 

The change in the power of (tt ) comes from a change 
in integration variable, required to match the dimension- 
less spatial input to F. Note also that V, which is not 
involved in the dynamics of the system, is left aside in 
the dimensional analysis. 

The resultant scaling function C, which describes all 
the unknowns with regard to this UETC, has no associ- 
ated absolute scale and is a function merely of two vari- 
ables. Further, since the quantities correlated are real 
(in real-space), then U* d (k) = U cd (—k) and statistical 
isotropy implies that the scaling functions are also real. 
Hence these functions are an efficient means of summa- 
rizing the data from the simulation as well as having the 
function of scale-extrapolation. 

It should be noted however, that the power of the 
UETC approach stems also from (i) strings decay on 
scales much smaller than the horizon and (ii) there can 
be no super-horizon correlations (since the strings form 
at the end of or after inflation). From (i), the string 
T^y (k, r) is unimportant for fcr 3> 1 and since the scal- 
ing functions involve a product of two such terms they 
therefore decay for high k\J tt' . In the opposite limit, 
(ii) implies that the scaling functions may be expanded 
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as simple power laws for low fcyr r' [45j]. Statistical 
isotropy then forbids odd terms in any such expansion 
(f/(k, r, t') = U(— k, t, t')) and the result is that 4 of the 
5 scaling functions required here are simply constant for 
low k\J t t 1 . In the final case, it is expected to vary as 
fc 2 rr', which can then be extracted to yield the same 
constancy for super-horizon scales. Furthermore, if the 
ratio of times t/t' is very large or very small, then the 
combination of (i) and (ii) implies that for any k^/r r', 
one of the factors of T M „ (k, r) is insignificant while one 
comes from the super-horizon regime. The scaling func- 
tions must hence decay for extreme time ratios and are 
most important for near-equal times. 

These simple considerations regarding the form of the 
scaling functions [53| mean that a field simulation, only 
able to provide a limited range in k\J t t' and t/t', can 
still provide the data required for accurate results over a 
large range of scales. 



B. Perturbation Equations 

As in the inflationary case, the perturbations evolve 
according to the linearized Einstein equations, covariant 
energy and momentum conservation and the Boltzmann 
equation [4||. The key change in the present case is the 
additional presence of the string energy-momentum ten- 
sor in the Einstein equations. In the linear regime ap- 
propriate for CMB calculations, the energy-momentum 
conservation equations may be separated into those in- 
volving the perturbations in the cosmic fluids and those 
involving the T^ v components, since products of these 
small quantities can be ignored. Hence, as stated in the 
Introduction, the defect energy-momentum tensor is sep- 
arately conserved. However, it also does not suffer from 
the gauge-dependence of the space-time metric and cos- 
mic fluid perturbations. This is as a result of it being 
identically zero in the homogeneous background relative 
to which the perturbations are measured. The metric 
and fluid perturbations are also described here in terms 
of gau ge- i nvariant quantities, via the formalism of Durrer 
et al. [iy, |45| to which the reader is referred for greater 
detail. 

It is convenient to decompose the perturbations into 3- 
scalars, 3- vectors and 3-tensors according to their trans- 
formations under 0(3) rotations. These three classes 
then evolve according to independent equations and 
are sourced by independent projections from the string 
energy-momentum tensor. Whereas vector perturbations 
decay in the standard inflationary scenario and so are 
not taken into account, they are continuously sourced in 
the defect case and must be considered. This is the sec- 
ond major difference between the defect and inflationary 
cases. 

As is common, the space-time metric is decomposed 
(in Fourier space) into scalar (S), vector (V) and tensor 



(T) parts as: 



.9oo 



(l + i s ) 



g 0i = a 2 (kiB s + Bj 



9ij = -a I 1 + SijH, 



k% kj q Oij 



(32) 
(33) 
(34) 



Hij 



In the scalar case, the gauge-invariant Bardeen potentials 
may then be formed as: 



3 a k 

k V a 



(35) 
(36) 



where S s = k^ 1 !!^ — B s . The Einstein equations for 
these potentials are then simplified if the two (of the four) 
scalar degrees of freedom in T M „ are projected out as: 



r.S 



■Ton 



qS qS rp 



3k m k n T n 



(37) 
(38) 



The relevant Einstein equations then appear as: 



fc 2 <l s = 4vrG(5|-3a 2 (p+p)$ s + ...) (39) 
k 2 (<l s + # s ) = 4ttG(5| + 5| + . . .) . (40) 



The ellipsis denote the additional presence on the right 
of the matter and photon perturbations, the notation for 
which shall not be defined here in a desire for brevity. 
Note that although these equations relate metric pertur- 
bations on the left to energy perturbations on the right, 
a Bardeen potential does appear on the right due the in- 
volvement of the metric tensor in raising and lowering 
the indices of tensor measures of the FRW background, 
namely the total background density p and pressure p. 
While there are actually four scalar degree of freedom in 
T^ v , that this tensor obeys covariant energy- momentum 
conservation provides two scalar equations via which the 
remaining two are specified. 

In the case of vector modes, it is useful to define: 



Toi kik m Tonn 



(41) 



which then obeys the vector constraint ki§Y= and so 
contains two of the four vector degrees of freedom in T^. 
This gives a simple form for the corresponding Einstein 
equation as: 



-k 2 £f= 16ttG 



sY 



(42) 



where l]7= 



k-'HY 



bY 



i . Momentum conservation then 
implies that the remaining two vector degrees of freedom, 
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which come from TV, , can be found from Si and its time 
derivative. 

Of the ten degrees of freedom in Ta V , the remaining 
two are 3-tensors and can be projected out of the space- 
space components in an analogous manner to ify . This 

T 

results in fiy, which sources perturbations via: 

Hij +2-H? +k 2 H l3 = 8ttG + ■ ■ •) ■ (43) 

There is no tensor equation that stems from energy- 
momentum conservation and hence this property cannot 
be used to apply further constraint in this case. 

The result is that the sourcing of perturbations by 
strings can be described in terms of the variables £§, 

S^, , Si and Sj , containing six degrees of freedom in to- 
tal. These appear in the Einstein equations, via which 
the space-time metric is determined, but do not other- 
wise influence the remainder of the standard CMBEASY 
evolution routine. That is with the exception that the 
continual sourcing of vector modes means that vector 
contributions must be additionally considered, and are 
evolved here using the Hu and White method [54j . 

C. UETC scaling functions 

Initially there are 10 independent components in the 
string energy-momentum tensor, between which there are 
potentially ±10 (10+1) = 55 UETC scaling functions to 
be considered. However as noted above, the number of 
scalar and vector degrees of freedom may be reduced 
via energy-momentum conservation such that there are 
two degrees of freedom for each scalar, vector and tensor 
class. Further, statistical isotropy enforces that UETCs 
between these classes are simply zero. 

For the scalar class, the three independent UETCs are 
EM: 

(§l(k, r) Si* (k, r'j) =-|Llcf 1 (kV^, r/r') (44) 

(St(k, r) S|*(k, r'j) = -0=lcf 2 (kV^, r/r') (45) 

(sf (k, r) £|*(k, r'j) =-?L±Cl 2 (kVr^, r/r'). (46) 

While the auto-correlations Cf x and C 22 are unchanged 
by the exchange of times r <-> r' (the scaling functions 
are real), the cross-correlation Cf 2 is not. Applying this 
transformation to Cf 2 gives the fourth member of this 
set, namely: 

(s|(k, r) S|*(k, r'j) = -iLl<5f 1 (W^ J r/r'), (47) 

and it is hence not independent: 

^(fcVTr 7 , r/r') = C?' (kVr^P, r'/r) (48) 



This will be calculated for the present work since it will 
always be the case that r > r' and hence the calculation 
of C 21 is required for the cross-correlation to be fully 
determined. 

In the vector case, the situation is slightly different in 
that the two independent degrees of freedom are spread 
across the three components of Si. It would clearly be 
inefficient to consider the correlations between all such 
components and hence the procedure chosen here is to 
project the two degrees of freedom by rewriting S Y in 
terms of an orthonormal basis (ki, e], e?): 

Sj= Hk^+elS^+efS^ 2 = e]S vl +e 2 S V2 . (49) 

The projection upon k is a scalar and is therefore zero, 
which leaves the projections upon e 1 and e 2 as the de- 
sired vector degrees of freedom. It is hence not necessary 
to explicitly calculate S^at all, but merely to apply this 
procedure directly to Ibi and discard the scalar projec- 
tion. 

Now, if all realizations are rotated 90° about a particu- 
lar k, then the two projections change as: 5 V1 — > ±S , 
S Y2 — > =fS vl . Hence statistical isotropy infers that their 
auto-correlations are equal: 

(s vi (k, r) S vl *(k, r'j) = (s V2 (k, r) S V2 *(k, r'j) (50) 

= eVr^^CVikVr'r 7 , r/r'), 

but that their cross-correlation is zero. In the present 
case, involving a finite number of realizations, the two 
auto-correlations are averaged. The above scaling func- 
tion definition matches that of Durrer et al. [U l45j |. 
although the basis projection approach here means that 
it is arrived at differently. 

It should be noted that the definition of CN differs from 
the general definition by a factor of k 2 rr' . This is relates 
to the fact that vector degrees of freedom can be written 
in terms of the curl of a vector field, and hence that 
Sy measures angular momentum. Therefore, the UETC 
must decay for large scales. As noted in the introduction 
to this section, at low k\Jrr' the scaling functions can be 
written as power laws expansions in k\J rr' and that they 
may contain only even terms. Hence in the vector case 
the scaling function may have a factor of k 2 rr' extracted 
so that all C tend to constants at low k\/rr' (4ll. |45|. 

Turning to the tensor modes, the source function Sij 
contains only two degrees of freedom and these may be 
projected out using a set of symmetric basis matrices: 

S7j= M^S 71 + M 2 S T2 . (51) 

Considering the vector basis above, it may be noted that 
hefef = and Sijefef = 8 AB . Hence two symmetric 
matrices which project out tensor modes are: 

Mij = ^=(ele 2 + e 2 e]) (52) 
M l = (53) 
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with the pre-factors required to make an orthonormal 
set. Consideration of a rotation by 45° is sufficient to 
show that the cross-correlations are zero while the auto- 
correlations are again equal: 

(£ T1 (k, r) S T1 *(k, r')) = (£ T2 (k, r) S T2 *(k, r')) (54) 
= 2^Ll(7T(fcV77 7 , t/t'). 

Note the factor of two, which is present to make G 1 have 
the same definition as the tensor scaling function used by 
Durrer et al. jMSl- 

D. Eigenvector decomposition 

Although the CMBEASY code produces the power 
spectra, the evolution of the perturbations is not per- 
formed using such quadratic quantities. Hence, although 
the UETC scaling functions defined above do contain all 
of the information required for CMB calculations, it is 
not immediately of the correct form for insertion into the 
source-enabled version of CMBEASY. However, suppose 
that a UETC set may be expanded as [10, H] : 

U bc (k-, t, t') = J2 A " {f «( fc ' T X( fc ' T ')- ( 55 ) 

n 

In the scalar case the indices b and c take on values 1 or 
2, with Ufi formed from Cf 1: but in the vector and ten- 
sor cases these indices are redundant. Now, the Green's 
function expression for the tensor component of a power 
spectrum (Eq. [23)) becomes: 

(x a (k,r)X* b (Krj) =J2 X ni:(k,r)ir(k,r), (56) 

n 

where: 

tf(fc,T)= J&T'g a {k,T,T)u n {k,T'). (57) 

As a result, the modified CMBEASY code can act upon 
the quasi-source functions u c n and so calculate one term 
in the sum, and hence eventually, the desired power spec- 
trum. 

Since scaling is broken at the radiation-matter transi- 
tion, then the UETC is given by a different scaling func- 
tion within each era. However, suppose that the ten- 
sor scaling functions calculated under both radiation or 
matter domination can be similarly decomposed. Numer- 
ically, C{kyTT' ', t/t') can only be represented at discrete 
values of its input parameters. Although it is not the 
form in which the simulations actually output data, a 
useful re-representation of the data is as C(kr, kr r ), for 
a set of discrete [kr]i and [kr']j with particular spacing. 
The scaling function is therefore represented as an M x M 
matrix Cij and the decomposition is then: 

Cy^A^e?. (58) 
n 



Being real and symmetric Cij is an Hermitian matrix 
whose eigenvectors form an orthonormal set. With that 
in mind, multiplication by c™ reveals that the supposed 
decomposition of Cij is nothing more than representing 
C^ in terms of its eigenvectors c™ and eigenvalues A m : 

CijC? = X m cT- (59) 

In the continuum, eigenvectors become eigenfunctions, 
but the discussion is little different, while numerically the 
decomposition can be straightforwardly achieved using, 
for example, the in-built functions of MATLAB [55| . 

If the scaling functions in the radiation and matter 
eras are similar in form (but not necessarily in mag- 
nitude), then the resulting eigenvectors will be similar. 
Hence a valid (but approximate) means of dealing with 
the radiation-matter transition is to write the integral /„ 
as [4l|: 

VKIn=f dr'g(k,T,r) -^L « x v/A^(fcT'), 

(60) 

where r denotes the radiation era and m the matter era. 
The function a r (T') is equal to one deep in the radiation 
era and zero deep in the matter era, with a parameter- 
izable transition between these limits as the dominant 
species changes (see [ll|). A similar approach is used to 
handle the transition from radiation to A domination, al- 
though no scaling functions are calculated in this era and 
the defects sources are merely allowed to decay (a x — > 0). 

While this discussion is only marginally changed for 
the vector and tensor contribution to the power spec- 
trum, for the scalar modes it is complicated by the in- 
volvement of multiple UETCs and the non-symmetric 
cross-correlation Cf 2 - In this case a 2M x 2M matrix 
is formed as: 

On = ( % 3 2 ) . (61) 

such that C^ = Cf 2 ([kT]i, [kr]j-M ) if i < M but 
j > M. This matrix is then symmetric on account that 
Cf 2 (kT, feT / )=(7| 1 (kr' , far) and these appear off-diagonal. 
From Eqs. 1551 and [551 the first half of the resulting eigen- 
vectors then represent the 5$, while the second half re- 
places Sty. 

Although CMBEASY does not actually involve such a 
Green's function approach, this eigenvector decomposi- 
tion allows the UETC data to be incorporated and the 
corresponding CMB power spectrum calculated as the 
sum of eigen-contributions. 

E. Collection of UETC data 

As already noted, the numerical simulations start from 
random initial conditions and therefore do not start in 
the scaling regime, or in fact even contain strings at very 
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early times. Hence there is an initial phase during which 
UETC data cannot be taken. Assuming that the system 
tends towards a scaling solution, eventually a time r sca Ung 
is passed at which the system scales to within the desired 
accuracy. The approach chosen here is to then store the 
six 3D arrays corresponding to the variables £§, S*|, S V1 , 
S V2 , S T1 and S T2 at that time. Then for each output 
time r > T sca ii ng , correlations are performed against these 
reference data. Angular averaging is performed (as well 
as the averaging of the two vector and two tensor UETCs) 
and the scaling functions Cf x , Cf 2 , Cfi , C22 > ^ an d 
output to disk as six ID arrays corresponding to the time 
ratio T/T sca ii ng . As previously noted, the output of (7fi 
is required since data for ratios less than unity are not 
calculated. 

Since the simulation is expected to be closest to scal- 
ing at later times, there is an argument for performing 
all unequal-time correlations against the final simulation 
time: U(k, r, r cn d). However, this means that the infor- 
mation required for correlation is not available until the 
simulation has been evolved to T enc i and the data from 
every prior output time must either be stored, or part of 
the simulation re-ran once the reference data has been 
acquired at r on( j. Given that the 3D arrays required for a 
512 3 lattice need 3 GB of storage, then storing them for 
each output time would be inefficient, while the second 
option requires the repetition of a significant proportion 
of the calculation. Additionally, it is the desire to run the 
simulations beyond their strict causal limit, in which case 
it is preferable to use such times merely for the less im- 
portant UETC data from the most extreme time ratios, 
rather than for the important equal-time data. 

IV. SIMULATION AND UETC RESULTS 
A. Simulations involved 

The results that are described here stem from 
15 000 CPU- hours of calculation on the UK National 
Cosmology Supercomputer [5l|, using 512 3 lattices and 
run across 64 CPUs. This enabled 5 realizations in both 
radiation and matter eras for s = 0.0, 0.2 and 0.3 and, ad- 
ditionally, s=l in the radiation era. At larger s values, 
the network formation is more challenging and the ini- 
tial period of non-scaling is generally larger. The largest 
value of s for which scaling could be attained quickly 
enough for the required UETC data to be taken from a 
matter era simulation was s = 0.3 and this is hence the 
closest value to the s = 1 goal. The s = 1 runs under radi- 
ation domination provide useful data on the non-artificial 
case, but scaling was not reached fast enough for all of 
the required UETC data to be collected and these are 
incorporated merely for reference. 

While of course this situation could be improved using 
larger (512 x M ) 3 lattices, which allow for longer causal 
run times, the processing time is very sensitive to such 
enlargements and varies as M . On the other hand, the 



actual increase in the available s values in the matter era 
would have been slight and in fact, from the results ob- 
tained, there is little to justify such a large computational 
outlay. 

B. Scaling in the string length 

It is obviously of great importance to test that the sys- 
tem approaches a scaling regime, and to a satisfactory 
level prior to the acquisition of UETC data. Further, 
the most important such test necessarily involves the ac- 
tual UETC itself, via which scaling is actually employed 
here. However as an auxiliary test, we additionally use 
the behaviour of the mean string length Lh in each hori- 
zon volume Vh- Given scaling, Lh grows in proportion 
to the horizon size r while Vh varies as t 3 , resulting in: 

C / « T . (62 ) 

V 

In the context of idealized-string simulations, this is 
expressed in terms of the Lorentz-invariant length and a 
simulation segment usually represents an element of this. 
In the present context the length must be derived from 
the fields and there is no unique means of proceeding. 
Vincent et al. [H| and Moore et al. [IIJ use a net wind- 
ing of the phase around the smallest closed loops repre- 
sented on the lattice, that is around lattice-plaquettes, 
in order to detect strings penetrating them. This enables 
the string paths to be traced and so the string length 
estimated. Here we employ a simpler approach and use 
the mean Lagrangian density in order to estimate the 
invariant length: 

^ = C ^ 

(since — C is the energy density for a static string, van- 
ishes for perturbative radiation, and is a 4-scalar). Al- 
though this is in many ways inferior to the above method, 
the results will be needed in the next sub-section and it 
will then be desirable that they directly involve the entire 
simulation volume. 

The results from 5 simulations in the radiation era with 
s = l and s~Q are presented in Fig. [3] The best-fit lines 
show a linear regime is reached of form: 

t cx (t-t c=0 ). (64) 

Similar results are found for the other simulations, in- 
cluding those in the matter era, with best-fit gradients 
as shown in Table HI It is encouraging that the slopes 
show no resolvable trend with s in cither era and that 
the s = case appears to be a good approximation to 
s = 1 under radiation domination (as was found previ- 
ously by Moore et al. [HI]). Further, in the s < 1 cases, 
which were ran with a final acausal period to extend the 
t/t' range, the trend continues to within the statistical 
uncertainties. 
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FIG. 3: Results for the Lagrangian measure of £ from 5 simu- 
lations in the radiation era with s = l (top) and s = (lower). 
The shaded regions show the 1 — a and 2 — a variations in the 
mean (error bars have not been used since correlations ex- 
tend across most of the plot and would tempt the reader into 
believing individual points were independent). The best-fit 
straight lines over the region 80</f> 1 < r < 128^o 1 (s = l) and 
6400 1 < r < 1280O 1 (s = 0) are also shown. Note that for 
the later case, this excludes the acausal final period. 



Simulation set 


£ versus 


r gradient 




radiation 


matter 


s = 1.0 


0.33 ±0.02 




s = 0.3 


0.299 ±0.014 


0.33 ± 0.03 


s = 0.2 


0.299 ±0.012 


0.304 ±0.012 


s = 0.0 


0.31 ±0.02 


0.304 ± 0.013 



TABLE I: The gradient of the linear region of the £ versus 
r plots from 5 realizations in each case. Note that results in 
this table do not involve the acausal over-run period. 

All runs show this offset scaling behaviour includ- 
ing smaller runs which employed a gauge-invariant (56| 
(but unsmoothed) version of the above phase- winding ap- 
proach to the string length. It is hence not a result of 
the string length measure employed but merely a con- 
sequence of the initial period of non-scaling. Choosing 
a different network formation history (that is a differ- 
ent r dependence of s at early times) gives no detectable 
change in the slopes, but has an effect on the r^ = o values. 
Hence, there is nothing fundamental about the offset and 
it is only the gradients that are important. Note that if 
the apparent linear regime was false, then it would be 
expected that the gradient too would be dependent upon 
the formation history (and therefore s as the early s(t) 
dependence must then be different). However, this is not 
the case. 

Further, if the seen trends are extrapolated to times 
of cosmological relevance, then the fractional disagree- 
ment between this and a direct proportionality becomes 
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FIG. 4: The raw equal-time scaling function Cfi(kr, kr) as 
averaged over 5 realizations for s = 0.3 in the radiation era. 
Results are plotted at roughly uniformly-spaced r values in 
the range 640^" 1 < r < 128^ 1 . The lower lines correspond 
to increasingly early times (dashed), with the 1 — a and 2 — <j 
uncertainties in the mean indicated for the latest time (solid). 
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FIG. 5: The equal-time scaling function Cfi(kr, kr) as in the 
previous figure, but with the time offset taken into account. 
For present plot, the mean offset across the 5 realizations is 
used to adjust results for each one, whereas the actual CMB 
calculations use independent offsets for each realization. 



insignificant and therefore these results are believed to 
be evidence for the required scaling behaviour. 



C. Scaling in the equal-time scaling functions 

A more complete test of the scaling hypothesis is af- 
forded by the UETC scaling functions, which should 
show no absolute temporal dependence, being a function 
merely of k\Jr t' and r/r'. This provides an analysis as 
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FIG. 6: The temporal dependence of the offset corrected 
CfiikTtkr) at four fixed values of kr: 17, 50, 110 and 220 
(which appear top — > bottom in the plot). Results are from 5 
realizations with s = 0.3 in the radiation era, using a common 
time offset for each realization. 

a function of scale and this is the test of direct relevance 
for the present CMB calculations. Since the functions are 
most important for equal (or near-equal) times, then a 
useful test is to compare the equal-time scaling functions 
calculated at different r. Typical such results are shown 
in Fig.g]for Cf ly during the period 640q 1 < r < 1280Q 1 
when £ varies linearly with r in the s < 1 cases. 

Unfortunately, a systematic variation is clearly visi- 
ble such that the overall normalization increases with r, 
while the peak becomes broader and shifts to the right. 
However, if all times in the calculation of C are replaced 
as: 

r = (r sim - 7£=o), (65) 

then the effect is a symmetric rescaling of the axes by 
a factor dependent upon the simulation time r s i m . This 
may be thought of as adjusting the results to estimate the 
scaling functions at late times, when the correction would 
have no effect. Once this process is performed, there is 
the desired time-constancy of the scaling functions, as is 
shown in Fig. [5] [69]. Similar results are seen for the other 
scaling functions and s values, but with slight exception 
of s=l. In this more ambitious case, good scaling is not 
seen for simulation times as early as r s - lm w 640^ 1 and 
reliable UETC data cannot be taken until later in the 
simulation. 

This behaviour may be further explored by plotting the 
C(fcr, kr) for a given kr as a function of time, as in Fig. [6] 
The results are qualitatively similar for the other scaling 
functions and s values, such that the important scales, 
close to the horizon, show scaling to a good approxima- 
tion before UETC data acquisition begins for the s < 1 
cases (r s j m = 64<^|~ 1 ). However, scales which are small 
relative to the horizon take longer to show convergence. 



Clearly this must be the case at some level in the present 
simulations since the string width is around one sixtieth 
of the horizon size, or larger, when UETC data is first 
collected. However, the use of such a dataset for CMB 
calculations should be accurate, considering the statis- 
tical uncertainties, given that the non-scaling effects at 
high fcr are unlikely to be larger than the statistical un- 
certainties at more important values. Hence, it would 
seem that the required UETC data is provided by the 
present s<l simulations with 512 3 lattices. 



D. Dependence of the ETCs upon s 

The dependence of such equal-time results upon s are 
shown in Fig. [7] for radiation era. This figure shows re- 
sults from the end of the causal period of simulation, by 
which time even the true s = 1 case shows good scal- 
ing. The results highlight that there is at most a very 
weak dependence upon s and that the typical difference 
between s — 1 and the artificial cases is approximately 
equal to the statistical uncertainties. Although an s = 1 
reference is not available in the matter era there are no 
obvious trends between the s < 1 cases for the scaling 
functions in that era, just as in the included figure. 



E. Unequal-time scaling function results 

Although the equal-time data does exhibit scaling, it 
is additionally required that the simulations are able to 
sample all of the important regions of the unequal-time 
scaling functions and, for example, the available r/r' is 
sufficient. This ratio is clearly constrained by the late 
onset of scaling, which is significantly later than the case 
of a non-linear cr-model, which has been previously used 
to represent global defects in field simulations [iol.l4ll.l47t. 
Here, we have additionally simulated global 0(4) textures 
via a very similar procedure to those previous works, in 
order to provide a useful reference with which to make 
comparisons against, but also as a check of the UETC 
and CMB calculation algorithms. In the texture case, 
the scaling regime may be comfortably studied in smaller 
256 3 simulations, with a maximum reliable time ratio of 
Rj 6. In comparison, the larger string simulations can 
provide a maximal ratio of 160/64 = 2.5 (s<l), which 
is then further reduced when the time offset is applied, 
becoming slightly less than 2 and dependent upon the 
precise case in question. 

The form of Cj^ under both radiation and matter dom- 
ination is shown for strings in Fig. [8l Results are shown 
for the s — 0.3 case but there are very similar for the 
other s values simulated: there is pronounced peak for 
t = t' and kr « 17, with a decline for unequal-times or 
for large k^/rr' values. Unfortunately, the late onset of 
scaling means that the super-horizon plateau is not well- 
covered by these simulations. It is hence shown here via 
an extrapolation at a constant level, determined from the 
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FIG. 7: The dependence of the C(kr, kr) upon s in the radiation era at T s i m — 1280 o 1 . The shaded regions show the 1 — a and 
2 — a uncertainties for the s = l case, while the lines indicate the s = 0.3 (solid), s — 0.2 (dashed) s = 0.1 (dot-dash) results. 




FIG. 8: The unequal-time scaling functions Cfi in the radiation era (left) and matter era (right). Results are from 5 realizations 
with s = 0.3. Each realization is given the average offset time such there is no interpolation involved which might otherwise 
smooth these plots. The data is then raw, albeit it has been extrapolated as a constant for low k\J t t' and the relative time 
symmetries have been used since the simulations output only for r > r' . Note there are large statistical uncertainties which 
cannot be shown here. 
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FIG. 9: The unequal time scaling functions from a global texture model, in order to provide a comparison case for the previous 
figure. The results show are from 10 realizations in the radiation era (left) and matter era (right), using a 256 3 lattice. Note 
the difference in the axis scales between this and the string case. 



mean of the function at each r/r' for the lowest k^/rr' 
values present. Further, it should be stressed that large 
statistical uncertainties exist which cannot be well-shown 
on these 3D plots. 

For comparison, the corresponding results from the 
texture model are shown in Fig. [9] It is immediately ap- 
parent that while the string case is dominated by a peak 
in both eras, this is not so in the texture case, which 
involves a merely small peak above the plateau in the ra- 
diation era (with a possible smaller and poorly resolved 
version under matter domination). However careful at- 
tention should be paid to the axes, which reveal that the 
string scaling functions are significantly larger than those 
of textures. Further, the string case shows significant 
contributions at relatively large k\J t t' values, whereas 
in the texture case decays for kr > 10 (t — t 1 ). Cos- 
mic strings are difficult structures to remove and while 
the highly curved super-horizon forms quickly straighten 
within the horizon, so that their length is very-much re- 
duced, strings do persist inside the horizon for some time. 
On the other hand, texture configurations may be con- 
tinuously smoothed away such that they are less signif- 
icant on smaller scales. This difference is highlighted 
by the sub-horizon peak in this string scaling function, 
which corresponds roughly to the inter-string separation. 
Clearly, these differences will manifest themselves in the 
resulting CMB predictions. 

The other unequal-time scaling functions from the 
s = 0.3 string simulations are shown in Fig. 1101 Each 
share an increase in magnitude and a shift to small scales 
relative to the texture case, although the corresponding 
plots are not shown due to brevity considerations. With 
the exception of the tensor functions, they are larger un- 
der radiation domination than in the matter era. How- 
ever, the overall forms do not change greatly between the 
two eras, which is desirable given the interpolation of the 
sources used in order to model the radiation-matter tran- 
sition. The differences are hence largely associated with 



the eigenvalues rather than changes in the forms of the 
eigenvectors. 

With regard to the question of r/r' coverage, all but 
and Cf 2 are clearly well-sampled, with significant 
decay at unequal times. It should be noted, however, 
that the cross-correlation function suffers from very large 
realization-to-realization differences and that the system- 
atic uncertainties arising from incomplete coverage of this 
function are likely to be insignificant relative to those of 
a statistical nature. This is also likely to be the case for 
the other functions and such systematic effects will be 
estimated as part of the CMB calculation. 



V. CMB POWER SPECTRA 

After the eigenvector decomposition of the above re- 
sults was performed, using a matrix size M = 512 and 
independent offsets for each realization, the modified ver- 
sion of CMBEASY was applied. This is the only point at 
which the cosmological parameters are involved and these 
were chosen to match the central values from non-CMB 
determinations: h = 0.72±0.08 ^bh 2 = 0.214±0.0020 
[H, n A = 0.75t ° 6 7 [HI; with additionally the inflation- 
motivated assumption of spatial flatness (and 0.1 used 
as the optical depth to the last-scattering surface). The 
string contribution to the temperature power spectrum 
is then given by the sum over eigen-contributions, each of 
which took of order 40 minutes of calculation time on a 
2.4GHz 64-bit AMD Opteron. This process is somewhat 
slower than the CMBEASY calculation of an inflation 
power spectrum since the vector mode must be addition- 
ally evolved, but also because the oscillating source func- 
tions necessitate a more careful integration of the equa- 
tions. The convergence of the eigen-sum in the s = 0.3 
case is shown in Fig. [Til with truncation after 128 terms 
giving convergence to within one percent. 

We estimate the statistical uncertainties in these power 
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FIG. 10: The remaining unequal-time correlation functions from s = 0.3 string simulations. For the cross correlation function 
C12 only the magnitude is shown and by virtue of Eq. [38] this is almost entirely anti-correlated. Note that its form can vary 
significantly from realization-to-realization and hence it suffers from very large statistical uncertainties. The vertical axes are 
constant within each scalar-vector-tensor class (including Fig. [8J so as to aid size comparisons. 
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FIG. 11: The convergence of the eigen-contribution sum for 
the temperature power spectrum due to cosmic strings. Re- 
sults are for s = 0.3 with (dashed) 1, 2, 4, 8, 16, 32, 64 and 
(solid) 128 terms included. 



FIG. 12: The CMB power spectrum contribution from cosmic 
strings simulated with s — 0.3 (solid), 0.2 (dashed) and 0.0 
(dot-dashed). The estimated 1 — a and 2 — a uncertainties (in 
the mean) are indicated in the s = 0.3 case by the shaded 
regions. 



spectra by repeating the CMB calculations using UETC 
results from individual realizations rather than first av- 
eraging the scaling functions. These are then shown in 
Fig. [T2]to be somewhat larger than the truncation error. 
Note however, that the statistical uncertainties are not 
directly related to cosmic variance since these results use 
the scaling approximation and do not stem from statisti- 
cal variations across a volume that actually corresponds 
to the observable universe. 

Additional systematic uncertainties stem from the lim- 
ited t/t' ratios achievable in the string simulations dur- 
ing the scaling regime. No extrapolation was made for 
the results presented here and the correlators were taken 
to be zero where data was unavailable. The correspond- 
ing systematic uncertainty can be explored by further 
zeroing all correlator data beyond a certain t/t' value, 
with the trend suggesting that the power spectra results 
are perhaps of order 10% too high as a result of the lim- 
ited r-range available, a result which is insensitive to the 
particular £ value. Finite-volume effects relating to the 
fact the system is only scaling approximately (particu- 
larly at earlier times in one or two realizations) can be 
explored by performing the £ fit over only a sub-set of 
the T s j m = 64 — ► 128 interval, ft is possible that the 
power spectra are underestimated by of order 10%, which 
again is not heavily dependent upon £. Other systematic 
sources of error include the matrix re-representation of 
the UETC data and associated numerical errors but a 
halving of the matrix size, for example, has a negligible 
effect on the results. The use of logarithmic rather than 
linear spacing in the fcr — fer plane slows the convergence 
of the eigen-contribution sum significantly, but upon con- 
vergence gives merely a slight shift on the high I of the 
peak and a change in the numerical results of 10% for 



I = 1000 but less than one percent for £ = 10, 100 or 
300. The modelling of the matter-radiation transition is 
a further source of uncertainty, for which an overestimate 
can be derived by the use of the matter scaling functions 
in the radiation era also. The contribution to large scales 
comes almost exclusively from the string sources in the 
matter era and the corresponding change in the power 
is negliblc for I ~ 10. The lower normalization of the 
matter era UETCs gives a reduction in power of 4% for 
£ = 100, which increases to 14% for £ = 1000, but these 
are clearly greater than the actual modelling errors. 

Fig. [12] additionally explores the effect of a change in 
s, which is of particular importance. It was initially an- 
ticipated that variations in s would provide a trend via 
which an extrapolation to the s — 1 result could be made. 
However, the £ and equal-time correlation results showed 
that any such variation was likely to be at most compa- 
rable to the statistical uncertainties, even for full range 
s= 0— >1 that was explored under radiation domination. 
The power spectrum results are consistent with this con- 
clusion and it appears that there is no basis upon which 
to perform any such extrapolation, with the variation 
from s = 0.0 to s = 0.2 being reversed for the jump to 
s = 0.3 and falling within the estimated statistical uncer- 
tainties. The position must therefore be taken that the 
present use of s < 1 simulations instills merely a system- 
atic uncertainty in the results, which is of a magnitude 
comparable to or smaller than the statistical uncertain- 
ties. For example, under a linear fit for the variation 
with s, the data of course allows for a zero gradient, but 
suggests an s — 1 extrapolation value that is 5% greater 
than that at s = 0.3 for £ = 10, although only 0.4% is 
due to the gradient, with the offset dominating. We refer 
the reader to the appendix for a more detailed discussion 
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FIG. 13: The decomposition of the CMB power spectrum 
(solid) into scalar (dashed), vector (dot-dashed) and tensor 
(solid-gray) modes. Results are shown for the s = 0.3 case. 



of the systematic uncertainties. 

The form of the power spectrum contribution is further 
studied for the s — 0.3 case in Fig. Q21 which includes 
a scalar-vector-tensor breakdown. This shows that the 
broad peak at 150 — 400 stems from both vector and 
scalar modes, which peak at l~ 180 and £«400 respec- 
tively. The power spectrum is hence dominated by vec- 
tor modes for all but the smallest scales, with the scalar, 
vector and tensor contributions having the approximate 
ratios 0.5 : 1.0 : 0.2 at ^=10. 

Although CMB data shows the actual power spectrum 
to peak at around £ ~ 200, with in fact a local mini- 
mum at £ « 400 [60l l6l| so that the string contribution 
must be sub-dominant, a useful comparison of forms is 
provided by setting the normalization of our results to 
match the WMAP 3-year data at I = 10, as in Fig. [H] 
Also included in this figure are the present global texture 
results, which match those from independent simulations 
0, EH EH and form a useful check of our algorithms, as 
well as providing a case for comparison. The slower de- 
cay of local strings within the horizon and the greater 
importance of their (scalar) UETCs in the radiation era 
yields a significant change in bias between high and low 
multipoles, with strings remaining important on smaller 
scales. The string contribution peaks when the data is 
close to a minimum and is also most precise. This may 
suggest that the fractional contribution from these string 
results would be more tightly constrained that those from 
textures [3j. However, it is also true that the overall form 
of the string contribution matches the data more closely, 
a fact that may over-turn such an argument. 

The normalization to the WMAP data at 
the conventional I — 10 (and the COBE mean 
temperature [62|) gives a value of 1-kG§\ as 
[2.04±0.06(stat.)±0.12(sys.)] x 10~ 6 , which is equal 
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FIG. 14: A comparison of the form of the s — 0.3 string 
contribution (solid) to the WMAP results [611 ]. enabled by an 
excessively large normalization to give a match to the data 
at £ = 10. Also shown are results from our global texture 
simulations (dashed). Note that the binned WMAP 3- year 
data are plotted which does not include output for £ = 10 
and that T is the mean CMB temperature . 

to G[i [631 ] - Given the sub-dominance of the string 
contribution, a value this large is actually ruled out by 
the data. Nevertheless, this number is useful for the 
comparison with results from alternative calculations, 
presented in the next section. Further, the result is not 
especially sensitive to the power spectrum normalization, 
varying merely as its square-root. Given the relative 
similarity in form with the texture result, then the 
power spectrum normalization needs to be reduced 
using a factor « 0.13 in order to fall in line with the 
first-year WMAP data incorporated in the likelihood 
analysis of Bevis et al. Q. However, this corresponds to 
a reduction of G/i using a factor of merely sa 0.4 and in 
the absence of a full likelihood analysis with the present 
string results (see [2l|), the WMAP-normalization value 
hence serves as guide-line upper-limit. 



VI. DISCUSSION AND CONCLUSIONS 

A. Comparison with previous string CMB results 

The key result presented is the form of the string con- 
tribution to the CMB temperature power spectrum, as 
shown in Fig. [14] with WMAP-normalization for com- 
parative purposes. Qualitatively, this form is similar 
to previous string determinations, including those from 
Nambu-Goto simulations [13] and the latest results from 
unconnected segment models (28[. There is agreement 
that the basic form of the power spectrum has a roughly 
constant slope at low multipoles, rising up to a sin- 
gle peak, with subsequent decay at small scales. How- 
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ever, the results from the present field simulations give 
a broader peak due to a greater relative contribution at 
large scales. While all three approaches yield results that 
peak at I ss 400, the idealized string methods each give 
75% of this peak value at I 100 whereas the present 
technique yields 75% of the peak value at a much lower 
multipole: I ~ 30. A comparison of the scalar-vector- 
tensor sub-contributions relative to the unconnected seg- 
ment results [HI , highlights that this dissimilarity stems 
largely from the difference in form of the vector contribu- 
tion. In the present case the vector component involves 
a very wide peak, quite different to that seen with the 
unconnected segment model. 

The £ = 10 WMAP-normalization of the power spec- 
trum result giving Gfi = (2.04 ± 0.13) x 10~ 6 
is larger than FRW Nambu-Goto simulation re- 
sults at large angular scales, for which COBE- 
normalization gave: Gu= (0.7 ± 0.2) x 10~ 6 [25j and 
Gfi = (imt° :w) x 10~ 6 [22] (with the dependence of this 
value upon the cosmology explored in the former). 
The UETC results of Contaldi et al. [H, involving 
Minkwoski space-time appears to be consistent with 
these, being G^ = 1.0x 10~ 6 (£ = 5). Finally this is 
true also for the latest unconnected segment result of 
Gn = l.l x 10~ 6 [H| (normalized to the total WMAP 
power but, by chance, an approximately equivalent nor- 
malization). 

It is difficult to be certain as to the cause of the differ- 
ences in form between the field evolution and idealized 
string results but possibilities include the treatment of 
the decay products or differences in velocity correlations. 
Although we consider the former to be more likely, it 
is interesting to note that a recent study of Nambu-Goto 
simulations in both Minkowski and FRW metrics [35[ has 
highlighted a difference with regard to velocity correla- 
tions and it is true that unconnected segment models 
cannot model such correlations. We hence advise some 
caution with regard to the interpretation of such results, 
which do not include the same level of physics as incor- 
porated in the present simulations. However, the extrap- 
olation involved here via scaling is over many orders of 
magnitude. Although it must be employed in all such cal- 
culations covering the full range of CMB scales and finds 
justification from the simulations themselves, its require- 
ment represents an inability to probe all scales involved 
to the desired degree. 



B. Comparison with global strings 

Global strings have previously been shown to give a 
contribution that is similar to that from global textures, 
albeit that the result peaks at slightly larger multipoles 
[i(j • Fig. [TJ] then indicates that global strings are inter- 
mediate in form between textures and the present local 
string results, but are actually somewhat closer to the 
former. Global strings have a less localized energy distri- 
bution than gauge strings and so experience significant 



long range forces [64| . This presumably causes the global 
case to show a more rapid decay within the horizon and 
their CMB contribution to have less significance at small 
scales, as seen in such calculations [40| . 



C. Future prospects 

It is not just the temperature power spectrum that is 
important. Significant constraints on cosmic string sce- 
narios might some day arise from the measures of the po- 
larization of CMB photons. Particularly, the so-called B- 
mode polarization spectrum provides an important win- 
dow on cosmic strings because inflation contributes to 
this only weakly. Scalar modes may contribute to the B- 
mode only via the gravitational lensing of the E-mode 
signal, with a second inflationary contribution arising 
from the sub-dominant tensor modes. It may hence be 
the case that the large vector contributions from cosmic 
strings enable their signature to be detected using data 
from future B-mode projects [65|, HH, [67| and this makes 
the difference in vector mode results especially interest- 
ing. Presently, B-mode results have only been published 
for unconnected segment models 28], with no direct in- 
put yet from Nambu-Goto simulations and therefore the 
situation is a little different in the polarization case. We 
will present polarization results from field simulations in 
a forthcoming publication [68j . 

The future also holds a great deal for the temperature 
power spectrum with, for example, the planned full-sky 
coverage at sub- WMAP scales from the Planck satellite 
[66| . Such data will help to restrict the inflationary con- 
tribution to the CMB and so more heavily constrain any 
sub-dominant component from cosmic strings. However, 
it should be noted that CMB perturbations from cosmic 
strings are non-Gaussian and are not statistically sum- 
marized using the power spectrum alone (differing in that 
respect from inflationary models under linear perturba- 
tion theory). Hence it is also important to consider pre- 
dictions for cosmic strings beyond the power spectrum, 
which our present UETC approach is unable to provide, 
but that may enable additional constraints. 
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APPENDIX A: SYSTEMATIC ERROR ANALYSIS 

In this appendix we detail the numerical experiments 
on which we base our error estimates, summarized in 
Sec. El The experiments are compared to our production 
runs, in which the parameters chosen are laid out in Table 

eh 

1. Conformal time offset tj = o 

The scaling form of the UETC functions (see Eqs. 
1441 - I55|) is extracted from the energy momentum cor- 
relators through multiplication by powers of the shifted 
conformal time r — r^ = o, as described around Eq. [64] 
Errors in the determination of the offset time t^ =0 will 
feed through into errors in the scaling functions. For 
the primary results, the offset time is found via a lin- 
ear fit of the string length parameter £(r) over the range 
64(^o 1 < T S i m < 1280Q . We estimate the errors by fit- 
ting over narrower ranges of conformal time, as shown in 
Table ITlTl These errors can also be viewed as an estimate 
of the errors from not reaching true scaling as t — > oo. 

2. Dynamic range 

Our data is taken over quite a limited dynamic range, 
defined as the maximum ratio R = t/t' , where r and r' 
are the two times in the UETC. Before the offset is taken 
into account, this ratio is 2.5, but is reduced to about 1.8 
with the offset (and is slightly different for each run). In 
order to estimate the errors from the limited dynamic 
range we further truncate the UETCs at smaller values 
of R, with results listed in Table ITVl We immediately see 
that the power spectrum increases with decreasing i? max , 
and that there is evidence that i? m ax ~ 1.8 overestimates 
the power. A logarithmic fit of the error against i? max 
indicates that, at least for I = 10 and £ = 100, the power 
spectrum is converging to a value about 10% lower than 
at -R max ~ 1-8. This is our estimate of the error due to 
the dynamic range. It is interesting that truncating the 
UETCs increases the power: this may well be due to the 
truncation producing "ringing" in Fourier space, which 
artificially sources extra perturbations. 

3. String width modification 

In order to deal with the shrinking of the string width 
in comoving units, we modify the equations of motion so 
that the strings grow in physical units. This growth is 
parametrized by s, the power of the scale factor by which 
the comoving width shrinks fEqs. [LA 1X5]) . As seen in ta- 
ble V and also in Fig. [7] there is no obvious trend to 
extrapolate and, as discussed in detail in Sec. V, the sys- 
tematic effect is likely to be comparable to the statistical 
uncertainties shown in the figure. 



Parameter Value 

Lattice size N 512 

Lattice spacing Aa: \(f>~ 1 ] 0.5 

Timestep At [fc x ] 0.1 

Scalar coupling Ao 2.0 

Gauge coupling eo 1.0 

Initial conformal time n [t^J 1 ] 1.0 

Final conformal time r e [(j)^ 1 ] 160 

Initial s -0.116 

Final s 0.3 

Time of change in s value [fl^g" 1 ] 32 

r sim range of £ fit [fa 1 ] 64 - 128 

Dynamic range _R max = (r/r') max « 1.8 

Eigenvector decomposition matrix size 512 

Eigenvector decomposition matrix kr spacing linear 

No. terms used in eigen-contribution sum 128 



TABLE II: Parameters used in production runs. The first 
part of the table lists the parameters of the Abelian Higgs 
model simulations, described in Sec. [TTJ with the second part 
listing parameters of the UETC method for calculating the 
CMB power spectrum as described in Sees. IIVI fVl 



r range 


I = 10 


I = 100 


£ = 300 


I = 1000 


64 - 128 














64 - 112 


-6.9 


-7.8 


-6.7 


-5.7 


64 - 96 


-8.7 


-9.4 


-8.6 


-9.7 



TABLE III: The percentage change in the temperature power 
spectrum when changing the conformal time range over which 
the string length parameter £ is fitted. 





£ = 10 


I = 100 i = 300 


I = 1000 


« l 


.8 











1.6 


7.0 


5.5 4.8 


8.1 


1.5 


13 


10 5.9 


11 


1.4 


22 


16 9.8 


14 


TABLE IV: The 


percentage change in the 


temperature 


power spectrum when changing the dynamic range parameter 


^ma X - 


= (r/r') max . 






S 




1 = 10 


I = 100 £ = 300 


£ = 1000 


0.3 













0.2 




13 


8.2 0.7 


-7.2 


0.0 




2.6 


-1.1 -2.0 


1.1 



TABLE V: The percentage change in the temperature power 
spectrum when changing the string width shrinkage parame- 
ter s. 





£ = 10 £ = 100 


£ = 300 


£ = 1000 


Matter & radiation 











Matter only 


0.03 -3.8 


-11 


-14 



TABLE VI: The percentage change in the temperature power 
spectrum when using matter era eigenvectors throughout the 
simulation, instead of interpolating to radiation era at r < 

7~cq • 
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No. eigenvectors 


^ = 10 


1 = 100 


£ = 300 


£ = 1000 


128 














150 


0.04 


0.09 


0.1 


0.2 



TABLE VII: The percentage change in the temperature power 
spectrum when changing the number of eigenvectors in the 
sum. 



Size of matrix 


I = 10 


£ = 100 


£ = 300 


I = 1000 


512 














256 


0.4 


0.2 


-0.4 


0.2 


128 


15 


-1.9 


-8.8 


-16 


64 


1.6 


-16 


-24 


-33 


TABLE VIII: 


The percentage change 


in the 


temperature 



power spectrum when changing eigenvector decomposition 
matrix size M. 



Type (No. e-vectors) 


£ = 10 


1= 100 


£ = 300 


I = 1000 


Linear (128) 














Log (128) 


-1.7 


-5.3 


-7.3 


-23 


Log (400) 


0.3 


0.8 


-0.8 


-11 



TABLE IX: The percentage change in the temperature power 
spectrum when changing eigenvector decomposition matrix 
sampling from linear to logarithmic spacing in kr. 



Source Value 

Conformal time offset +9% 

Dynamic range —10% 

String width modification ±5% 

Matter-radiation ±11% 

Eigenvector sum < 1% 

Eigenvector decomposition matrix size < 1% 

Eigenvector decomposition matrix kr spacing < 1% 

Total in quadrature ±16% 



TABLE X: Estimates of systematic errors in the temperature 
power spectrum at £ = 300 due to the sources discussed above. 

4. Matter-radiation transition 

It is not obvious in the UETC method what to do when 
the expansion rate changes. We interpolate between ra- 
diation era and matter era eigenfunctions. In order to 
estimate the error associated with this procedure we cal- 
culated the power spectrum with matter era eigenfunc- 



tions only, which had the effect of decreasing the power 
by up to 14% at i = 1000, but was negligible at £ = 10 
(see table EH- 



5. Eigen-contribution sum convergence 

Fig. [TT] shows how the power spectrum changes as in- 
creasing numbers of eigenvectors of the UETC matrices 
are included. To estimate how close we are to conver- 
gence we increased the number of eigenvectors from 128 
to 150, finding the negligible changes given in Table IVlTl 



6. Eigenvector decomposition matrix size and 
spacing 

The UETCs are recorded at discrete values of kr and 
r = t/t' , but must be diagonalized as a matrix with 
rows and columns labeled by kr and kr' . The number 
and spacing of the kr (and kr') values used to construct 
this matrix are not fundamental and must be chosen to 
give reasonably accurate results while requiring a mini- 
mal number of eigen-contributions for the convergence of 
the power spectrum. It seems that the size of the ma- 
trix barely affects the power spectrum above a value of 
256, as can be seen from Table IVIIII We tried linear 
and logarithmic sampling, finding that there was negli- 
gible difference providing enough eigenvalues were taken 
in the sum, as shown in Table iLXl In view of the smaller 
number of eigenvalues linear sampling is preferred, and 
we conclude that negligible errors are associated with the 
sampling. 



7. Summary 

The important range for future fitting to WMAP data 
will be I k, 300, for which we summarize the errors in 
Table IXl Note that the conformal time offset error gives 
an underestimate while the dynamic range error means 
that our results are likely to be an overestimate. We ob- 
tain final estimates for the upward and downward errors 
by combining these with the other errors in quadrature. 
When applied to the WMAP normalization value of G/i 
we find Gfi = [2.04 ± 0.06(stat.) ± 0.12(sys.)] x 10~ 6 , re- 
calling that the temperature power spectrum depends on 
the square of G/i. 
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tions would not have provided data for the earliest and 
latest times. Given the already short range of times over 
which scaling may be examined and a desire to keep the 
results presented as raw as possible, independent offsets 



were not used until the eigenvector decomposition phase, 
when the interpolation was already needed in order to 
transfer the UETC data onto the matrix grid. 



